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By implementing the nonlinear Poisson-Boltzmann theory in a cell model, we theoretically investigate the influence of poly elec¬ 
trolye gel permeability on ion densities and pH deviations inside the cavities of ionic microcapsules. Our calculations show that 
variations in permeability of a charged capsule shell cause a redistribution of ion densities within the capsule, which ultimately 
affects the pH deviation and Donnan potential induced by the electric field of the shell. We find that semipermeable capsules 
can induce larger pH deviations inside their cavities that can permeable capsules. Furthermore, with increasing capsule charge, 
the influence of permeability on pH deviations progressively increases. Our theory, while providing a self-consistent method for 
modeling the influence of permeability on fundamental properties of ionic microgels, makes predictions of practical significance 
for the design of microcapsules loaded with fluorescent dyes, which can serve as biosensors for diagnostic purposes. 


1 Introduction 

Polyelectrolyte (PE) microcapsules - ionic (charged) colloidal 
particles with hollow cavities 1 2 - have attracted great atten¬ 
tion in the past decade due to their novel fundamental prop¬ 
erties and their potential applications as biosensors to moni¬ 
tor local ion concentrations (such as pH) in cellular environ¬ 
ments. 3 14 Charged PE shells, which may encapsulate pH- 
sensitive fluorescent dyes, generate electric fields that can 
cause deviations in local ion distributions and limit practi¬ 
cal applications. Previous experiments 1516 and theories 1718 
demonstrated that the pH and ion densities near a charged 
flat surface can deviate from their bulk values. For exam¬ 
ple, Bostrom et al . 15 demonstrated that ion and pH gradients 
emerge near biological flat membranes. Zhang et alJ^ ex¬ 
perimentally confirmed that local ion concentrations and sen¬ 
sor read-outs can be attributed to surface charges. Janata^ 17 18 
showed theoretically that pH shifts measured by optical sen¬ 
sors depend on bulk-surface interactions. The charged shell of 
a spherical microcapsule can induce variations in ion densities 
near the strongly curved surface of the capsule. Furthermore, 
because of the asymmetry between the inside and outside en¬ 
vironments, the measured ion concentrations in microcapsule 
cavities can deviate greatly from those in bulk. Understanding 
such ion deviations is significant for biomedical applications 
of ionic microcapsules as biosensors, e.g., to_avoid misdiag¬ 
nosis of diseases, such as early-stage cancer] 19 21 Ion distri¬ 
butions can also affect the release properties of pressurized 
capsules for drug delivery. 22 

Ion density deviations inside of microcapsule cavities have 


a Department of Physics, North Dakota State University, Fargo, ND 58108- 
6050, USA. E-mail: alan.denton@ndsu.edu 

b Current address: Institutfiir Theoretische Physik, Georg-August Universitat, 
37077 Gottingen, Germany 


been studied previously by experiments,®! 6 theories,!£5|27| 
and simulations. EH For example, Sukhorukov et al 2 ^and Gao 
et al . 24 demonstrated that a pH difference between the inside 
and outside of PE capsules emerges when the capsules are per¬ 
meable to small ions, but exclude poly(styrenesulfonate) ions 
of a particular molecular weight. Another approach, based on 
entrapping polyanions within PE microcapsules, also reported 
a redistribution of H + ions across a semipermeable microcap¬ 
sule wall. 25 These ion density differences were attributed to 
the capsule semipermeability and explained by a macroscopic 
theoretical model of Donnan equilibrium. 25 27 Vinogradova et 
al. 29 31 applied Poisson-Boltzmann (PB) theory to study the 
osmotic pressure acting on semipermeable shells in poly ion 
solutions. Stukan et al . 28 used molecular dynamics simulation 
to study ion distributions near nanocapsules that are permeable 
to solvent and counterions, but impermeable to poly electrolyte 
coils, which were modeled as soft colloids. Semipermeability 
to salt ions has also been linked to osmotic shock of rigid pro¬ 
tein shells, such as viral capsidsJ 32 

Recently, we demonstrated, using non-linear PB theory, that 
even fully permeable microcapsules can induce ion density de¬ 
viations inside of charged capsule cavities. 33 These deviations 
depend on the degree of dissociation of the PE making up the 
shells, rather than on the permeability of the capsule wall. In 
practical applications of PE microcapsules in cellular environ¬ 
ments, however, the capsule walls usually exclude poly ions, 
such as charged DNA or amino acids. In this case, the ion 
density deviations induced by the microcapsule should also 
depend on the permeability of the charged wall. The influ¬ 
ence of permeability on deviations of local ion densities and 
pH induced by ionic microcapsules is still poorly understood. 

In this paper, by modeling a PE microcapsule as a uniformly 
charged shell that is permeable to some ionic species (counte¬ 
rions and salt ions), but only semipermeable to another species 
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(polyions), we extend our previous work to analyze the impact 
of permeability of charged capsule shells on ion density devia¬ 
tions in aqueous solutions. By employing nonlinear PB theory, 
we systematically calculate deviations of local ion densities 
inside microcapsule cavities induced by the permeability of 
capsule shells. Unlike macroscopic models of Donnan equi¬ 
librium for neutral capsules, 27 here the redistribution of ion 
densities depends not only on the capsule wall’s permeability, 
but also on its degree of dissociation (charge density). Our 
results demonstrate that capsule permeability can also signifi¬ 
cantly influence the properties of ionic microcapsules and their 
performance as biosensors. 

The remainder of the paper is organized as follows. In 
Sec. [2] we define our model of semipermeable ionic micro¬ 
capsules and describe our implementation of PB theory. In 
Sec. [3] we present numerical results from our calculations for 
the influence of capsule permeability on ion densities, Donnan 
potentials, and pH distributions. Finally, in Sec. [4| we sum¬ 
marize our findings and emphasize implications for practical 
applications. 


2 Models and Methods 

2.1 Uniform-Shell Model of Ionic Microcapsules 

We consider a bulk dispersion of ionic capsules, microions, 
and polyions dispersed in water. Within the primitive model 
of poly electrolytes,El the solvent is idealized as a dielectric 
continuum with uniform dielectric constant £. We model the 
capsules as spherical shells of inner radius a , outer radius b , 
and valence Z. In a polar solvent, a capsule becomes charged 
when counterions dissociate from the PE chains that form the 
cross-linked network (hydrogel) making up its shell. Simi¬ 
larly, the polyions become charged by dissociation. In line 
with many experiments, we presume the capsules to be nega¬ 
tively charged. Approximating the distribution of ionized sites 
in the hydrogel as uniform within the volume of the shell, we 
represent the number density of fixed charge by a simple radial 
profile, 


n f {r) 


( 0 , 


4;r(Z? 3 — a?) 

,0, 


r > b , 

a<r<b , 
0 < r < a , 


( 1 ) 


(a) 



(b) 



Fig. 1 (a) Model of a spherical microcapsule of inner shell radius a, 
outer shell radius b, and valence Z in water. The capsule is permeable 
to water and small ions (smaller spheres labeled as + and —), but only 
partially permeable to polyions (larger, shaded spheres), (b) Bulk 
dispersion of microcapsules in Donnan equilibrium with a reservoir. 


ther sign, depending on the species and the pH level of the 
solution. For simplicity, we model microions and polyions as 
point charges. For reasons made clear below, we furthermore 
consider only monovalent ions. 

The capsules are permeable to water and microions, but 
only partially permeable to polyions. To quantify the degree 
of penetration of the capsule shell by polyions, we introduce 
a permeability factor a, which ranges from 0 (no penetration) 
to 1 (complete penetration). Within the cavity (r < a ), we 
take the dielectric constant to be the same as in bulk, while in 
the shell, we assume a lower dielectric constant (£ s h e ii < £), 
as suggested by experiments on PE microgels in water] 35 ® 
In Donnan equilibrium, the capsules are confined to a fixed 
volume, while the solvent, microions, and polyions can freely 
exchange with a reservoir, with fixed number densities of salt 
ion pairs, no, and of polyions, n p 0 [see Fig.[Jb)]. 


2.2 Poisson-Boltzmann Theory of Bulk Dispersions 


We model a bulk dispersion of ionic capsules via Poisson- 
Boltzmann theory. 37 In its density-functional formula¬ 
tion, 38 39 PB theory focuses on the grand potential functional, 
£l[n±(r),n p (r)], which is a unique functional of the num¬ 
ber densities of positive and negative microions, n±( r), and 
of polyions, n p ( r). We arbitrarily assume positively charged 
polyions of valence z p . In reduced form, the PB approxima¬ 
tion for the grand potential functional can be expressed as 


p£l[n±(r),n p (r)] = Jdr in + ( r) 


In 


ft+( r ) 

no 


-1 


where r is radial distance from the center of the shell 
[Fig.[lja)]. Although neglecting charge discreteness and ion- 
specific effects, this coarse-grained model is consistent with 
the primitive model and valid on length scales longer than the 
typical spacing between ionized groups on the PE backbones. 

The microions comprise dissociated counterions and salt 
ions (e.g., Na + and Cl - ). The polyion charge can be of ei- 
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+ il dr [n+ (r )-n- (r) + z p n p (r )-n f (r)] \jf(r) , (2) 

where the first integral accounts for the ideal-gas free energy 
and the second integral for the electrostatic potential energy. 
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Here, /3 = 1 /(k B T) and i/r(r) = /Je0(r) is a dimensionless 
form of the electrostatic potential 0(r), 


e(j)(r) = Jdr'[n + (r)-n-(r)+z p n p (r)-nf(r)]v(\r-r '\), 

( 3 ) 

generated by Coulomb pair interactions, v(r) = e 2 /(er ), be¬ 
tween elementary charges e on mobile microions dispersed 
throughout the system and fixed ions localized within the 
shells. The neglect of correlations among ions inherent in the 
mean-field approximation for £1 is valid for weakly-correlated 
monovalent ions, but questionable for more strongly corre¬ 
lated multivalent microions J 29 31 For this reason, we restrict 
our considerations to monovalent ions. 

In thermodynamic equilibrium, the grand potential func¬ 
tional is a minimum with respect to the microion and polyion 
densities. Minimizing Eq. ^ with respect to n±( r) and n p ( r) 
yields Boltzmann distributions for the equilibrium densities: 


and 


n±( r) = noe^ 





-ZpVW 


a n-ne z P^( r ) 


r>b , 

0 < r < b , 


( 4 ) 

( 5 ) 


where the factor a dictates the permeability of the capsule to 
poly ions, ranging from completely impermeable for a = 0 to 
completely permeable for a = 1. Substituting Eqs. ([4]) and ^ 
into the Poisson equation, 


V 2 <Kr) = - T p(r), 


( 6 ) 


where 


p(r)=e[n+(r)-n-(r)+z p n p (r)-n f (r)\ (7) 


is the local charge density, yields the Poisson-Boltzmann 
equation: 

vV(r) = K$swh\if(r)+47tkB[nf(r)-z p rip(r)] . (8) 

Here, X B = e 2 j ( £k B T ) is the Bjerrum length (A# = 0.714 nm 
in water at temperature T ~ 293, £ = 80); K*o = \/%7iX B nQ is 
the inverse Debye screening length in a reference reservoir 
of pure salt solution (absent polyions); and K p = yJ\nX B n p r) 
can be interpreted as an inverse Debye screening length of the 
polyions. By solving Eq. ^ with appropriate boundary con¬ 
ditions, we obtain the microion and polyion density profiles. 


2.3 Poisson-Boltzmann Cell Model 

In general, solving the PB equation [Eq. © with boundary 
conditions matching an arbitrary arrangement of capsules is 
numerically quite challenging, although feasible in some sys¬ 
tems, such as charged colloids. 40 42 For computational effi¬ 
ciency, we instead use a cell model, 43 44 whose boundary con¬ 
ditions are relatively simple. The cell model is justified by the 
wide disparity in size and charge between capsules and mi¬ 
croions and by our focus on solutions of relatively low salt 
concentration E5E3 

For spherical capsules, the cell model represents a bulk 
dispersion by a spherical cell - centered on a single cap¬ 
sule - of radius R determined by the capsule volume frac¬ 
tion rj = ( b/R ) 3 [see Fig. [lja )]. Along with the capsule, 
the cell contains counterions and salt ions, which may freely 
penetrate the capsule, and polyions, whose penetration of the 
capsule is limited by the shell’s permeability. The condi¬ 
tion of electroneutrality of the cell relates the ion numbers: 
Z = N+-N-+ZpN p . 

In the spherical cell model, the PB equation simplifies to 


y/'(r) + -y/(r) = < 


Kq sinhy^(r) — z p K^e ZpV(r) , 

^ (yo sinh y(r) - azp^-^^ + ) . 

Kq sinh \jf(r) - az p K 2 e~ z P ¥ ^ , 


b<r<R , 

a< r<b , 
0 < r < a , 


( 9 ) 


where r is the radial distance from the center of the cell and 
X = £sheii/£ < 1 is the ratio of the dielectric constant in the 
capsule shell to that in the bulk solvent. 

Boundary conditions on Eq. © impose continuity of the 
electrostatic potential at the inner and outer boundaries of the 


capsule shell, 

V^in(^) = Vshell(^) 5 V^shell (^) = tyout {b) , (10) 

vanishing of the electric field at the center of the cell and on 


This journal is ©The Royal Society of Chemistry [year] 


Journal Name, 2010, [volf 1-[9] | 3 











the cell boundary, 


K( o) = o, i|/„„ t (fi)=o, (ii) 

as required by spherical symmetry and electroneutrality, and 
continuity of the electric displacement on the inner and outer 
shell boundaries, 

Vm(a) = ZV4di( fl ) , xV s h e ii( b ) = Yout( b ) , (12) 

where the solutions in the three regions are labelled as 1//i n (r) 
(0 < r < a ), V'sheiiM (a < r < b), and i/wM (b < r < R). 

By numerically solving the PB equation [Eq. 0]. along 
with the boundary conditions [Eqs. (^Q|-(^2|)], in the three ra¬ 
dial regions (inside the cavity, in the shell, and outside the 
capsule), we calculate the equilibrium microion and polyion 
density distributions within the spherical cell. 

3 Results and Discussion 

3.1 Distribution of Electrostatic Potential and Field 

To illustrate our theory, we present numerical results for 
system parameters typical of experiments on microcapsules. 
Specifically, we consider negatively charged microcapsule 
shells of inner radius a = 50 nm, outer radius b = 75 nm, and 
valence Z = 500 dispersed in water at room temperature (A# = 
0.714 nm), in Donnan equilibrium with a salt reservoir of con¬ 
centration no = 0.1 mM, at volume fraction Tj =0.1, corre¬ 
sponding to a cell radius R = ri~ l / 3 b ~ 3.23 a. This concen¬ 
tration is sufficiently dilute to ensure independence of the ion 
distributions within neighboring cavities. We set the dielectric 
constant ratio between the microcapsule shell and the solution 
at X = 0.5, which is consistent with measured dielectric con¬ 
stants in hydrated ionic Poly(N-isopropylacrylamide) (PNI- 
PAM) microgels ranging from 63 at 15°C to 17 at 40°C . 35 E3 
In order to justify the mean-field PB approach, we consider 
here only weakly correlated monovalent polyions (z p = 1). Fi¬ 
nally, we choose the polyion concentration, n p o = OAno, to 
give a polyion screening length K~ l somewhat longer than the 
screening length in a reference reservoir of pure salt so¬ 
lution. 

Figure [2ja) shows the distribution of electrostatic potential 
within the cell. Starting from the center, the potential de¬ 
creases with increasing radial distance, reaching a minimum 
within the shell. The depth of this minimum increases with 
decreasing capsule permeability, reaching a value of ^1 k^T 
(in potential energy) for a = 0.1. The kinetic barrier to ther¬ 
mal diffusion of ions across the shell is thus sufficiently low 
to ensure equilibrium ion distributions. Outside the shell, the 
electrostatic potential increases toward the cell edge. 

The corresponding electric field is shown in Fig.|2|b). Start¬ 
ing from zero at the center, the field decreases as r increases 



Fig. 2 (a) Distribution of reduced electrostatic potential \\f(r) vs. ra¬ 
dial distance r from center of a microcapsule of inner shell radius 
a = 50 nm, outer radius b = 75 nm, valence Z = 500, and dielectric 
constant ratio % = 0.5 in an aqueous solution of microcapsule vol¬ 
ume fraction 17 =0.1, reservoir salt concentration no = 0.1 mM, and 
polyion concentration n p = 0.04 mM. Results are shown for shell 
permeability a = 0.1, 0.5, and 0.9, (b) Electric field E(r) vs. r for 
permeability a = 0.1. Inset shows distribution of displacement field 
D(r) = eE(r). 
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inside the cavity, but then rises sharply within the shell. The 
discontinuities at the shell boundaries (r = a and 1.5 a) orig¬ 
inate from the difference in dielectric constant between the 
shell and the solution. The displacement field, rather than 
the electric field, is continuous at the boundaries [see inset 
to Fig. [2jb)]. Outside the shell, the field decreases as r in¬ 
creases, approaching zero at the edge of the cell to match the 
outer boundary condition. As a consistency check on our im¬ 
plementation of the theoretical model, we calculate the total 
charge number N tot of microions and polyions in the cell, 


f R 

N tot =4nj^ drr 2 [n + (r)-n-(r)+Z p n p (r)} , (13) 

and confirm global electroneutrality (N tot = Z) over a range of 
permeabilities. 

3.2 Donnan Potentials 

Chemical equilibrium between the inside and outside environ¬ 
ments of ionic microcapsules can result in a significant dif¬ 
ference in electrostatic potential between the two sides of the 
shell, as revealed in Fig.|2ja). We define this difference as the 
shell Donnan potential'. 

Yd (shell) = y(b) - <//(«) , (14) 

which is analogous to the Donnan potential at the surface of a 
bulk polyelectrolyte gel.® The variation of the shell Donnan 
potential with capsule permeability a is shown in Fig. [3ja). 
With increasing permeability, polyions increasingly penetrate 
the capsule shell, leading to a gradual decrease of ty/p(shell). 
For comparison, Fig. [3jb) shows the Donnan potential at the 
cell edge, y/z)(cell) = !//■(/?), as a function of capsule perme¬ 
ability. Interestingly, i//£>(cell) is relatively insensitive to vari¬ 
ation of a, implying that capsule permeability also does not 
significantly influence the bulk osmotic pressure of the disper¬ 
sion. On the other hand, the permeability is expected, to in¬ 
fluence ion distributions inside and outside the capsules, with 
potential impact on practical applications of ionic microcap¬ 
sules, such as in diagnosis of diseases. In the remainder of this 
section, we examine in detail ion density distributions induced 
by shell permeability, as well as corresponding pH deviations 
inside microcapsule cavities. 




0.0 0.2 


0.4 0.6 0.8 

a 


Fig. 3 Donnan potential \f/o vs. capsule permeability a at the edge 
of (a) the capsule shell (r = a) and (b) the entire cell (r = R). Other 
system parameters are as in Fig. [2] 


3.3 Influence of Permeability on Ion Densities 

Next, we investigate the influence of capsule permeability on 
the ion density distributions within the cell. In practice, per¬ 
meability may increase upon swelling, triggered by a change 
of solution pH or temperature. Figure |4ja) shows the distri¬ 
bution of polyion density over a range of permeabilities. As 
the capsule becomes more permeable (a increases), poly ions 
can more easily penetrate the capsule shell, thus increasing 
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the polyion density within the cavity. This redistribution of 
polyions results from entropy-driven thermal diffusion, which 
can overcome the electrostatic potential energy barrier [see 
Fig.0a)]. 

The rising concentration of polyions inside the capsule 
upon increasing permeability leads to an expulsion of positive 
microions and a corresponding decrease in n+(r) for r < b, 
as shown in Fig.[4jb). Nevertheless, there is a net increase in 
the total concentration of positive ions inside the capsule. To 
maintain electroneutrality of the system, the density of nega¬ 
tive microions also increases, as shown in Fig. (Jc). The re¬ 
distributions of ion densities induced by changes in capsule 
permeability should be mirrored by corresponding changes in 
the local concentrations of H + and OH - ions, leading to devi¬ 
ations in local pH from bulk solution values. 

3.4 Influence of Permeability on pH in Microcapsules 

Next, we investigate the influence of microcapsule permeabil¬ 
ity on pH deviations inside the capsule cavities and explore the 
potential impact on applications to pH sensors in cellular envi¬ 
ronments. The polyions in our model are positively charged, 
corresponding to the real scenario of amino acids dispersed 
in an alkaline cellular environment. Under such conditions, 
the pH within the cell is determined by the concentration of 
hydroxyl (OH - ) ions. Assuming the concentration of OH - 
ions to be proportional to that of all negative microions, the 
local pH deviation induced by the charged shell can be ap¬ 
proximated by 

ApH(r) =log[n-(r)/n-(R)] . (15) 

The average pH deviation inside the capsule cavities is then 

(ApH) = ^ [“ drr^ApU^r) . (16) 

a* Jo 

For a capsule of valence Z = 500, Fig.[5]shows (a) the local pH 
deviation profile and (b) the average pH deviation induced by 
the charged shell over a range of capsule permeabilities. With 
increasing a, the local and average pH deviations inside the 
cavity decrease in magnitude, as a result of variations in ion 
concentrations discussed in Sec. |3.3[ Thus, in alkaline envi¬ 
ronments, pH deviations inside cavities of ionic microcapsules 
induced by the charged shells are suppressed by increasing 
permeability to positive polyions. In other words, semiperme- 
able charged capsules can induce larger pH deviations inside 
their cavities that can permeable capsules. 

In a previous paper, 33 we demonstrated that the valence of 
ionic microcapsules dominates pH deviations inside the cavi¬ 
ties. The larger the capsule valence, the greater the pH devi¬ 
ations. Here we find that the permeability of the capsule can 
also influence pH deviations. Finally, we investigate how the 
deviations induced by permeability vary with capsule valence. 



r/a 



Fig. 4 Density distributions of (a) polyions, (b) positive microions, 
and (c) negative microions vs. radial distance r from center of micro¬ 
capsule for permeability a = 0.1, 0.5, and 0.9. Other system param¬ 
eters are as in Fig. [2] 
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Fig. 5 (a) Deviation of local pH from bulk value vs. radial distance r 
from center of microcapsule [Eq. (JT5J] for valence Z = 500 and per¬ 
meability a = 0.1, 0.5, and 0.9. (b) Deviation of average pH inside 
cavity (r < a) vs. permeability [Eq. |I6|]. Other system parameters 
are as in Fig. [2] 
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Fig. 6 (a) Average deviation of pH from bulk value inside micro¬ 
capsule cavity vs. valence Z for permeability a = 0 and 1. (b) 
Permeability-induced average pH deviation vs. valence [Eq. GZh 
Other system parameters are as in Fig. [2] 
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Figure[6ja) illustrates the variation with capsule valence of the 
average pH deviation inside the cavity for the two extremes of 
permeability, a = 0 and 1. We find that, with increasing va¬ 
lence, the average pH deviation increases in magnitude. For 
capsules that are impermeable to polyions (a = 0), however, 
the increase is more rapid than for fully permeable capsules 
(a = 1). To quantify this effect, we define 


7= (ApH) a=0 — (ApH) a=1 , (17) 


which measures the average pH deviation inside the capsule 
cavity induced solely by permeability of the capsule shell. 
Figure |6jb) shows the permeability-induced average pH de¬ 
viation inside the cavities as a function of capsule valence, il¬ 
lustrating that, as valence increases, the average pH deviation 
associated with capsule permeability increases in magnitude. 
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4 Conclusions 

In summary, by applying a cell model implementation of the 
nonlinear Poisson-Boltzmann theory to solutions of ionic mi¬ 
crocapsules that are semipermeable to polyions, we have an¬ 
alyzed the influence of capsule permeability on ion densities 
and pH deviations inside the capsule cavities. In real biolog¬ 
ical systems, variations in microcapsule permeability could 
result from variations in porosity of the polyelectrolyte net¬ 
works making up the capsule shells associated with the size 
distribution of polyions (e.g., amino acids) or from network 
swelling/de-swelling induced by changes in temperature, pH, 
or salt concentration. 

Our results show that, upon varying capsule permeability, 
the ion densities redistribute so as to fulfill the competing re¬ 
quirements of minimum free energy and global electroneutral¬ 
ity. Ultimately, increasing permeability suppresses deviations 
in microion density and pH induced by the charged shells. 
These findings have potential relevance for the design of mi¬ 
crocapsules that encapsulate fluorescent dyes to serve as ionic 
biosensors for diagnostic purposes. Although we have fo¬ 
cused, in this study, on solutions containing only monovalent 
ions, in order to justify our use of the mean-field PB theory, 
our approach could be extended to multivalent ion solutions 
by incorporating ion correlations into the PB theory! 47 ! or by 
performing molecular simulations in the cell model. 
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